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ABSTRACT 

Direct numerical simulation of turbulent homogeneous shear flow is performed in order to 
clarify compressibility effects on the turbulence growth in the flow. The two Mach numbers 
relevant to homogeneous shear flow are the turbulent Mach number M t and the gradient Mach 
number M g . Two series of simulations are performed where the initial values of M g and M t 
are increased separately. The growth rate of turbulent kinetic energy is observed to decrease 
in both series of simulations. This ‘stabilizing’ effect of compressibility on the turbulent energy 
growth rate is observed to be substantially larger in the DNS series where the initial value of M g 
is changed. A systematic comparison of the different DNS cases shows that the compressibility 
effect of reduced turbulent energy growth rate is primarily due to the reduced level of turbulence 
production and not due to explicit dilatational effects. The reduced turbulence production is not 
a mean density effect since the mean density remains constant in compressible homogeneous 
shear flow. The stabilizing effect of compressibility on the turbulence growth is observed to 
increase with the gradient Mach number M g in the homogeneous shear flow DNS. Estimates 
of M g for the mixing layer and the boundary layer are obtained. These estimates show that 
the parameter M g becomes much larger in the high-speed mixing layer relative to the high- 
speed boundary layer even though the mean flow Mach numbers are the same in the two flows. 
Therefore, the inhibition of turbulent energy production and consequent ‘stabilizing’ effect of 
compressibility on the turbulence (over and above that due to the mean density variation) is 
expected to be larger in the mixing layer relative to the boundary layer in agreement with 
experimental observations. 


Supported by National Aeronautics and Space Administration under NASA Contract No. NAS1-19480 while 
the author was in residence at the Institute for Computer Applications in Science and Engineering (ICASE), 
NASA Langley Research Center, Hampton, VA 23681-0001. 


1 



1. Introduction 


The inhibited growth of the shear layer thickness and the turbulent stresses in the high- 
speed mixing layer is a well-known phenomenon (Bradshaw 1977; Kline, Cantwell and Lilley 
1982; Papamoschou and Roshko 1988). For example, the thickness of a plane mixing layer 
at convective Mach number M c =1.5 grows at a rate which is only a third of the nominal 
incompressible value. However, the reasons for the strong ‘stabilizing’ effect of compressibility 
in the turbulent mixing layer remain unclear. One explanation (Sandham and Reynolds 1989; 
Morris, Giridharan and Lilley 1990) draws an analogy with the result of linear analysis that the 
growth rate of small disturbances decreases when the convective Mach number M c increases. 
Here, the convective Mach number M c is the ratio of the mean velocity difference (U\ — U 2 ) 
across the mixing layer to the sum (ci + c^) of the sound speeds in the two streams. Linear 
stability analyses (Sandham and Reynolds 1989; Ragab and Wu 1989; Jackson and Grosch 
1990) show that the growth rate of the most amplified disturbance decreases as a function of 
M c . Jackson and Grosch (1990) have shown that the linear analysis results of the maximal 
growth rate for a wide range of values for the free-stream temperature and density fall on 
essentially a single curve which is a function of M c . The similarity between the variation of 
disturbance growth rate as a function of M c , and the variation of the experimentally observed 
growth rate of shear layer thickness as a function of M c has prompted the argument that linear 
stability analysis explains the compressibility effect of reduced mixing. However, it is unclear 
how and why does simple linear theory apply to the fully turbulent mixing layer. 

Another explanation of the reduced mixing of the high-speed mixing layer has been ad- 
vanced by Zeman (1990) and Sarkar et. al. (1991a), who suggest that the dissipative effect of 
dilatational velocity fluctuations (such fluctuations have V ■ u / 0) becomes progressively more 
important when the turbulent Mach number M t increases, reduces the turbulent energy, and 
thereby decreases turbulent mixing. The turbulent Mach number M t = u/c is the ratio of the 
r.m.s. (root mean square) velocity fluctuation u to the speed of sound c. Turbulence models 
developed by Zeman (1990) and Sarkar et. al. (1991a) which parametrize the dilatational corre- 
lations as functions of M t have been able to capture the decreased growth rate of the high-speed 
mixing layer. However, the hypothesis that dilatational effects lead to the stabilizing effect of 
compressibility needs direct validation. 

Compressibility effects in the high-speed boundary layer are different from those in the high- 
speed mixing layer. When the free-stream Mach number M ^ of the boundary layer increases, 
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there is a moderate decrease in skin friction and thickness growth rate which is directly due to 
the reduction of density from its free-stream value and the consequent reduction in momentum 
transport by the Reynolds shear stress. In the case of the compressible mixing layer, an increase 
in the convective Mach number M c leads to a dramatic reduction in the growth rate of the shear 
layer thickness, a reduction which is far larger than can be explained by the density variation. 
The discrepancy between compressibility effects in free shear layers and wall boundary layers 
is not understood. A theoretical explanation for the reduced growth rate of the high-speed 
mixing layer should account for this crucial difference between turbulent free shear flows and 
wall-bounded flows in the high-speed regime. 

Due to the limitations of current computer hardware, three-dimensional DNS of the com- 
pressible mixing layer or the compressible boundary layer have not been performed at sufficiently 
high Reynolds number at which the turbulence would be fully developed. It is relatively easier 
to simulate homogeneous turbulent flows at realistic turbulence Reynolds numbers, and further- 
more perform a parametric DNS study that compares results from a few different cases. In the 
present paper, simulations of homogeneous shear flow, a flow with constant mean shear rate, 
are described. Homogeneous shear flow retains the feature of sustenance of turbulent kinetic en- 
ergy by mean shear common to a variety of turbulent shear flows. The homogeneous shear flow 
problem has been studied previously with DNS by Feiereisen et. al. (1982), and more recently 
by Sarkar, Erlebacher and Hussaini (1991b), and Blaisdell, Mansour and Reynolds (1993). Al- 
though the two latter DNS investigations did find that increased levels of compressibility lead to 
decreased levels of turbulence just as observed in experiments on the high-speed mixing layer, a 
systematic comparison of different cases to identify possible causes for this ‘stabilizing’ effect of 
compressibility was not performed. Furthermore, the studies neglected to discriminate between 
the gradient Mach number M g and the turbulent Mach number M t . The gradient Mach number 
is M g = Sl/c, where S = dU / dy is the mean shear rate, / a representative integral length scale 
of the turbulence in the direction of shear, and c the speed of sound. The parameter M g can 
be viewed as the ratio of an acoustic time l/c for a large eddy to the mean flow time scale 1 / S . 
Durbin and Zeman (1992), in their RDT (rapid distortion theory) analysis of homogeneous 
compressed turbulence, recognized that a parameter similar to the gradient Mach number that 
is defined by Am = Dl / c with D = |V • U\ is one of the relevant parameters; however, these 
authors considered only the limit Am << 1. Later, Jacquin, Cambon and Blin (1993), and 
Cambon, Coleman and Mansour (1994) showed, again for homogeneous compressed turbulence, 
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that compressible RDT results are sensitive to Am. It should be noted that the RDT restriction 
of high mean distortion rate relative to the lare-scale turbulence distortion rate does not apply 
to the present simulations. 

Since the gradient Mach number M g introduced in the present work is based on the mean 
velocity held like the convective Mach number M c , it is important to compare and contrast 
these two Mach numbers. There is of course a similarity, because M g = Sl/c can be viewed 
as a ratio of the mean velocity difference SI across a ‘large-scale eddy’ to the speed of sound, 
just as M c is related to the ratio of the mean velocity difference across the two streams in a 
mixing layer to the speed of sound. However, there are essential differences. The parameter M g 
is a held quantity which varies across an inhomogeneous shear layer unlike M c . Moreover, M g 
differs between shear hows such as the mixing layer, boundary layer, and the wall jet since the 
variation of mean velocity held and turbulence length scale is different in these hows even if the 
overall velocity difference across the hows is the same. In fact, we will show later in this paper 
that the value of M g is much larger for the supersonic mixing layer relative to the supersonic 
boundary layer even if the mean how Mach number has the same value in the two hows. 

In our earlier DNS studies (Sarkar, Erlebacher and Hussaini 1991b) of compressible homo- 
geneous shear how, the speed of sound was varied keeping all other parameters hxed which led 
to the simultaneous change of initial gradient Mach number M g q and initial turbulent Mach 
number M t o- Similarly, M g q and M t o were simultaneously increased (or decreased) in the in- 
dependent DNS investigation of Blaisdell, Mansour and Reynolds (1993). The present study 
performs two series of new simulations that comprise six different cases to address the following 
questions: (1) How do the gradient Mach number and turbulent Mach number individually 
effect the how evolution? and (2) What are the probable reasons for any compressibility effects 
on the how evolution due to the variation of these Mach numbers? In series A of the simula- 
tions, M g o is varied keeping M t o constant and vice-versa in series B. It should be noted that 
M g and M t are different parameters because their ratio M g /M t = SI /u is not a constant, that 
is, SI /u can depend on other factors such as the how geometry, compressibility, wall effects and 
initial/boundary conditions. Thus among incompressible hows, the log region in the bound- 
ary layer, equilibrium homogeneous shear how, and the near-wall region in the boundary layer 
have progressively larger values of Sl/u and correspondingly different turbulence structure and 
growth rates. For example, Lee, Kim and Moin (1990) have shown that incompressible 
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homogeneous shear flow with initial Sl/u that is sufficiently larger than the equilibrium value 
exhibits low- and high-speed streaks just as in near-wall turbulence. 


2. Preliminary Analysis 

The compressible Navier-Stokes equations provide the mathematical model for the problem. 
After performing the Reynolds decomposition into mean and fluctuating parts of the velocity, a 
system of equations can be obtained for the evolution of the instantaneous density p*, velocity 
fluctuation v,;*', and the instantaneous pressure p*. The superscript ★ is used to denote a 
dimensional variable, the overbar is used to denote an averaged variable, and the superscript 
/ is used to denote fluctuations in a variable with respect to its average. The overtilde ~ 
will be used later for denoting mass-weighted or Favre average. The following equations are a 
specialization of the general equations to the homogeneous shear flow problem where the mean 
velocity is T^ = Sx £, where S is a constant. 


dtP * + U tp* + (p*u*')j — 0 (1) 

d t (p*ui*’) + + (/ W')j = ~P*i ~ Sp*u 2 *'Sii + (2) 

dtp k + U*jPj + K V),j = -(7 - 1 )p k uf tj + (7 - 1)k2^ + (7 - 1)$* (3) 

p* = p *BT* (4) 


where 4>* = + pS 2 is the dissipation function, p the molecular viscosity, R the gas 

constant, 7 the ratio of specific heats, and k the thermal conductivity. The fluid properties p, 
R and k are taken to be constant for simplifying the problem. The fluctuating viscous stress is 


= K u tj 


+ u %i) 



( 5 ) 


We wish to determine the variables that parametrize compressibility effects on the turbu- 
lence evolution. One method is to determine the conditions under which the acoustic time scale 
Iq/cq becomes important for the behavior of a given turbulence velocity fluctuation uq with 
length scale Iq. For this purpose, we non-dimensionalize equations (??)-(??) with Iq/cq as the 
time scale, the initial mean density po as the density scale, the initial mean temperature To as 
the temperature scale, uq the fluctuating velocity scale, To the mean velocity scale, and Iq the 
length scale. The instantaneous pressure is nondimensionalized by PqUq. The exception is the 
first term on the right hand side of Eq. (??) which is explicitly split into a mean pressure and 
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fluctuating pressure, and the mean pressure is non-dimensionalized by the initial mean pressure 
Po so that P = P* /P 0 = 0(1). The nondimensionalization adopted here uses a hybrid scaling 
because incompressible scaling po M o i s use( i f° r pressure variation while acoustic scaling l/c 0 is 
used for the time variation. 

The nondimensional equations become 


dtp + M 0 U ip : i + Mto{ pui'\i — 0 


(6) 


dt(pu i) + MqU j(pu/)j + M t o(puj'ui)j = -Mtop'i - Mgopu^Sn + (7) 

dtp+MoU^ + Mtoiu/p), = - - ( 7 -l )M to p'u' jjj + -^-( p / p ) tjj + 


jM t0 


Pr 0 Re 0 


Ren 


l^toP 

p 


(8) 

(9) 


The non-dimensional parameters unrelated to the speed of sound that appear in Eqs.(??)- 
(??) are the Reynolds number Re o = uolo/vo and the Prandtl number Pro = pC p jn. Since 
the mean convection term that contains U i can be removed in the special case of homogeneous 
turbulence by transforming the equations to a frame moving with the mean velocity, the mean 
Mach number Mq does not effect the evolution of the flow variables in that moving reference 
frame. Thus, the only two ‘compressibility’ parameters in the equations for homogeneous shear 
flow are the turbulent Mach number M to = uq/cq and the gradient Mach number M g o = SIq/cq. 
For a general flow, the gradient Mach number can be defined by M g q = Slo/co where S = 
\J Ui,j Ui,j , and Iq an appropriate integral length scale. In the case of shear flow, Iq is chosen to 
be the integral length scale of the stream wise fluctuating velocity in the shearing direction X 2 - 

Equation (??) suggests that the fluid momentum may change significantly on the acoustic 
time scale if M t o = 0(1) or M g o = 0(1). Similarly, Eqs. (5) and (7) show that, if M t o = 0(1) 
the variation of fluid density and pressure on the acoustic time scale can be significant. Thus, 
it is clear that compressibility effects on the flow evolution increase when either M g q or M t o 
increases. In the present DNS study, the quantities M g q and M t o are considered to be the 
compressibility parameters and varied between different cases. 

We note that the above discussion of compressibility effects on the flow evolution is relevant 
only after any initial transients in the flow have subsided. The role of initial conditions in the 
evolution of compressible isotropic turbulence has been previously investigated by Erlebacher, 
Hussaini, Kreiss and Sarkar (1990) with a DNS of the full system and with a low M t o asymptotic 
analysis where the nonlinear and viscous terms were dropped. The reduced equations in the 
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analysis constitute a hyperbolic system whose eigenvalues show strong asymmetry essentially 
due to the l/M t o scaling of the right hand side of Eq. (??) versus the M t o scaling of the right 
hand side of Eq. (??). We found that, due to the asymmetry in the eigenvalues, dilatational 
velocity fluctuations or pressure fluctuations can grow rapidly on the acoustic time scale from 
their respective initial values under some circumstances, and furthermore classified the condi- 
tions which lead to the initial transients. The initial conditions in this study have been chosen 
so as to remove the influence of such initial transients. 

The turbulent kinetic energy K = u\u\j 2, where the Favre-average u\u\ = pui'ui' /~p is 
an estimate of the energy associated with the velocity fluctuations. For the homogeneous 
turbulence considered here, K is computed as a volume average over the computational domain. 
The turbulent kinetic energy is a function of time and the non-dimensional parameter A(/) = 
(1 / S K)(dK / dt) describes the instantaneous value of its temporal growth rate. It may be 
expected that physical effects that decrease A in homogeneous shear flow will inhibit turbulent 
energy levels and thereby inhibit shear layer growth in inhomogeneous turbulent shear flows as 
well. In fact, it is shown below through an analogy that the coefficient C$, in the commonly 
accepted expression for the streamwise growth rate of a mixing layer between two streams with 
uniform speeds U\ and U 2 

dS_ ffi - U 2 
dx S U\ + U 2 ’ 

is proportional to A. 

Consider a self-similar, plane mixing layer between two streams of velocity U\ and U 2 - The 
kinetic energy K in such a mixing layer is given by 

K = Coilh - U 2 ff(r ] ) (11) 


where Co is a constant, the similarity variable r] = y/S(x ), and S(x) is the shear layer thickness. 
Let K* be the kinetic energy integrated across the transverse direction, that is, 


K* 



/ +00 

fid) dr] 

-OO 


From Eq. (??), it can be shown that 


1 dK* _ 1 d6_ 
K* dx S dx 


( 12 ) 


(13) 
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We now consider homogeneous shear flow, U(y) = Sy with S a constant. The temporal 
growth rate of the turbulent kinetic energy K(t) is described by the nondimensional variable A 
defined by 


1 dK 
SK dt 


(14) 


where S is the constant shear rate. Physical experiments (Tavoularis and Karnik 1989) and 
DNS (Rogers, Moin and Reynolds 1987) support the notion that in the case of incompressible, 
homogeneous shear flow, the turbulence eventually grows exponentially, that is, A(/) — ► A^ 
(a constant) for large t. Although this picture of turbulence is consistent with available data, 
there is no rigorous proof that the turbulence grows exponentially for all time. In any case, A 
is a useful nondimensional estimate of the instantaneous growth rate of K. Assume that the 
temporal evolution of the turbulent kinetic energy can be transformed into a spatial evolution 
of K* by the transformation x = U c t , where U c is a characteristic convection velocity given by 


u c = C 2 (U! + U 2 ) 


(15) 


where C 2 is a constant, U\ and U 2 are the velocities at the top and bottom of the shear layer, 
respectively. Then Eq. (??) can be transformed into a relation for the equivalent spatial growth 
rate 

1 dK* 

K* dx 


where A^ is the constant, large-time value of the kinetic energy growth rate. 

Now we draw an analogy between the homogeneous shear layer and the mixing layer, and 
equate the right hand side of Eq.(??) and Eq. (??) to obtain 

dS _ Aq 0 lh- U 2 
dx C 2 U\ + U 2 

According to the analogy, the growth rate of the mixing layer thickness is proportional to the 
asymptotic growth rate of turbulent kinetic energy in the corresponding homogeneous shear 
flow. Furthermore, Eq. (??) implies that a constant value of A^ in the homogeneous shear 
flow problem is consistent with the linear growth rate, dS jdx oc (U\ — U 2 )/(Ui + U 2 ) (Schlichting 
1979), seen in fully developed, constant density, mixing layers if the coefficient C 2 in Eq.(??) 
for the convection velocity Uc is a constant independent of the free-stream velocities. 


AooA 


U c 

Aoo (Ui ~ U 2 ) 
SC 2 (U 1 + U 2 ) 


(16) 
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In summary, this section establishes that the gradient Mach number and the turbulent Mach 
number should be considered as important compressibility parameters in homogeneous shear 
flow. Furthermore, it is shown that the behavior of the normalized temporal growth rate of 
the turbulent kinetic energy, A = (1 / SK)(dK/dt), as a function of the gradient and turbulent 
Mach numbers is important because of its connection to shear layer growth in more general 
turbulent shear flows. 

3. DNS method 

The algorithm for the DNS of compressible homogeneous shear flow is essentially identical 
to the spectral collocation, third-order Runge-Kutta algorithm which was used in our previous 
simulations. Details can be found in our previous work (Sarkar, Erlebacher and Hussaini 1991b). 
All the simulations were performed on a 128 3 spatial grid in a cube of side 27T. Each simulation 
had the same initial data which corresponded to uniform density p = 1, solenoidal velocity 
fluctuations, pressure fluctuations calculated from the usual incompressible Poisson equation, 
and temperature obtained from the equation of state of an ideal gas. Such a choice minimizes the 
buildup of compressibility effects due to initial transients from incompatible initial conditions. 
The initial velocity is a Gaussian random field whose energy spectrum E(k) oc k 4 exp( — 2 k 2 /k 2 m ) 
with k m = 18. The flow evolution is governed by the following parameters which need to be 
prescribed: the shear rate S , the viscosity //, the thermal diffusivity k/ poC p , the initial mean 
speed of sound Co, the initial r.m.s. (root mean square) velocity uq , and the initial turbulence 
length scale Iq. The corresponding nondimensional parameters are M g o, M t o, Re\ o, and Pro 
which are computed as follows: M g o = Slo/co , where S is the constant mean shear rate, Iq 
the initial integral length scale of u in the transverse shearing direction and cq the initial mean 
speed of sound; M t o = uq/co , where uq = \^2Kq and Kq is the initial turbulent kinetic energy; 
Re\o = MoAo/zz, where the Taylor microscale is Ao = uo/ujo, and uq = yEJicJi is the r.m.s. 
vorticity; and Pro = pC p / k. 

Two series of simulations labeled A and B, respectively, were performed. In series A, the 
gradient Mach number M g o is progressively increased in cases A1 to A4 by changing the shear 
rate S , all other non-dimensional parameters remaining the same. In series B, the turbulent 
Mach number M t o is progressively increased in cases B1 to B3, all other non-dimensional 
parameters remaining the same. Both S and c were appropriately changed in series B to so 
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as to vary M t o while M g o simultaneously remains constant. Tables 1 and 2 give the non- 
dimensional parameters applicable to the simulations performed here. Note that Case A1 is 


Case 

MgO 

Afto 

Re\ o 

Pro 

(SK/e) 0 

A1 

0.22 

0.40 

14 

0.7 

1.8 

A2 

0.44 

0.40 

14 

0.7 

3.6 

A3 

0.66 

0.40 

14 

0.7 

5.4 

A4 

1.32 

0.40 

14 

0.7 

10.8 

Table 1 

: Parameters for series A of the DNS 

Case 

MgO 

Afto 

Re\ o 

Pro 

( SK/e) 0 

B1 

0.22 

0.13 

14 

0.7 

5.4 

B2 

0.22 

0.20 

14 

0.7 

3.6 

B3 

0.22 

0.40 

14 

0.7 

1.8 


Table 2: Parameters for series B of the DNS 


the same case as Case B3 so that there are a total of six different simulations. 

The range of variation of M g o and M t o is chosen to cover a reasonable portion of parameter 
space with case B1 as an example with low gradient and turbulent Mach numbers, and case A4 
as an example with high gradient and turbulent Mach numbers. The constant value M t o = 0.4 
in series A is chosen to be somewhat high so that (SI /u ) o = M g o/M t o does not become so large 
as to cause the turbulence evolution to be dominated by linear effects (rapid distortion theory). 
The initial value (SI /u ) o varies between 0.55 and 3.3 in series A. The Prandtl number is chosen 
to be Pro = 0.7 since the flow of air is considered here. The large scales of turbulence acquire 
energy due to the mean shear while the energy of the small scales grows due to the increase in 
turbulence Reynolds number and consequent increase in nonlinear energy transfer. Thus, the 
integral length scale / increases and the Kolmogorov length scale r] decreases with time. The 
simulations were terminated when either r] became small enough to cause insufficient resolution 
of the small scales by the grid, or when / became large enough for the periodic boundary 
conditions to affect the growth of the large scales. The termination time of the simulations 
varied between St = 15 and St = 20. The turbulent Reynolds number increases with time and 
its value at the end of the simulations varied between Re\ = 35 and Re\ = 42 which indicates 
that the velocity held becomes turbulent, albeit with moderate Reynolds number. Since the 
turbulence is spatially homogeneous, turbulence statistics and correlations at a given time are 
obtained by averaging over the 128 3 points in the computational domain. 
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4. DNS results on the growth of turbulent energy 

The compressibility effect on the evolution of turbulent kinetic energy K = u'-u'-/2 is dis- 
cussed here. The non-dimensional parameter A = (1 /SK)(dK/dt) describes the temporal 
growth rate of the turbulent kinetic energy. It was shown in section 2 that the behavior of A is 
of fundamental interest. 

The equation for turbulent kinetic energy K is now analyzed to identify the factors that 
influence the growth rate. The exact equation for the evolution of K in homogeneous shear flow 
is, 

C ^- = V - e s - e c + p’d’ /p (18) 

where V = — Su^u^ is the production, e s = z/adad the solenoidal dissipation rate, e c = (4/3)z7d /2 
the compressible dissipation rate, and p'd 1 the pressure-dilatation. Here lo\ is the fluctuating 
vorticity, and d' is the fluctuation in dilatation (V • u). It is useful to work with an evolution 
equation for the growth rate A = (1 / S K)(dK / dt), an equation which is obtained from Eq. (??) 
to be, 


A 


— 2&12 - 


^ S 

Jk 


(e c - p'd' /p) 

SK 


-26i 2 (1 


£ s 

v 


(e c - p'd'/p) 

SK 


(19) 


where 612 = u' x u' 2 I2K is the Reynolds shear stress anisotropy. The last term in Eq. (??) repre- 
sents the explicit influence of compressibility through dilatational effects. In the incompressible 
case (Speziale 1991), it is generally accepted that homogeneous shear flow turbulence evolves 
to a state with constant values of 612 and e s /V , and therefore, from Eq. (??), A also approaches 
a constant equilibrium value A^. It is clear from Eq. (??) that a change in the equilibrium 
value Aqq in homogeneous shear flow with increasing compressibility can be due to any of the 
following: (1) a change in the magnitude of 612 which implies a more (or less) ‘efficient’ pro- 
duction mechanism, (2) a change in e s /V , the relative dissipation of the turbulence, and (3) a 
change in (e c — p'd')/SK, which represents dilatational effects. 


4.1 Results from Series A 

Fig. 1 shows the values of A at integral values of St in the simulation. It is clear that there 
is a systematic decrease of the asymptotic (large St) value A^ of the turbulence growth rate 
when M g o increases in cases A1-A4. For example, case A4 with M g 0 = 1.32 has A^ ~ 0.03, 
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which is more than a factor of 3 smaller than the corresponding value of A^ ~ 0.11 in case A1 
with M g0 = 0.22. 

The terms on the right hand side of Eq. (??) are now evaluated and compared in order to 
identify possible causes for the reduced growth rate of K. Fig. 2 shows the evolution of the 
nondimensional production (V / SK = —2612). By St = 20, the turbulence eventually evolves to 
approximately constant values of Reynolds shear stress anisotropy in each case. Furthermore, 
the long-time values of 612 show a systematic decrease with the corresponding value of M g 0. 
The magnitude of — 2612 in case A4 is smaller by a factor of 3 relative to that in case Al. Thus, 
the efficiency of turbulent production by the mean shear is strongly inhibited when the gradient 
Mach number increases. 

The evolution of e s /V is shown in Fig. 3a. After an initial transient, there is only a slight 
difference between the different cases which is far too small to cause the large reduction in A. 
For example, e s /V at St = 20 varies between 0.5 and 0.6 for cases A1-A4. The 20% increase 
in in the relative dissipation e s /V in case A4 relative to case Al is a small effect compared to 
the factor of 3 decrease in the nondimensional turbulent production ( — 2612). Fig. 3b shows 
the evolution of the normalized dissipation e s /SK. The long-time values of e s /SK show a 
systematic decrease from case Al to A4. Since e s /SK = —2(e s /V)bi2, such a decrease is 
consistent with compressibility decreasing the value of — 612 (in Fig. 2) while having very little 
effect on e s /V (in Fig. 3a). 

Fig. 3c shows the behavior of the dilatational terms. Although small, the normalized dilata- 
tional term in Fig. 3c is not negligible; its value at the end of the simulation ranges between 
10-30% of the growth rate A of the corresponding case. Due to a similar observation in our ear- 
lier simulations (Sarkar, Erlebacher, and Hussaini 1991b), we attributed significant importance 
to dilatational terms. However, the systematic comparison of the different DNS cases conducted 
in the present study shows that the dilatational terms do not lead to the compressibility effect of 
decreased growth rate of turbulent kinetic energy. This can be seen by comparing the different 
curves in Fig. 3c. The difference in the normalized dilational terms between the various cases 
is small and also of the wrong sign to explain the difference in the values of A^ between these 
cases. For instance, the values for the dilatational term (e c — p'd')/ SK are smaher for St > 15 
in case A4 relative to case Al, implying that dilatational effects would increase the growth rate 
in case A4 relative to case Al contrary to the observed decrease in A. 
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Eq. (??) for A can be rewritten as 


A — —2612(1 — X t 


( 20 ) 


where the influence of all terms other than the production term is lumped into X t which is 
given by, 




e s + e c - p'd'/p 
V 


( 21 ) 


Fig. 4 shows that, after the initial transient, X t becomes approximately constant and shows 
very little difference for St > 15 between the different cases. Thus, the large reduction in the 
value of A is almost wholly due to the large reduction in the magnitude of Reynolds shear stress 
anisotropy 612. 

For homogeneous flow, the Helmholtz decomposition gives a unique decomposition of the 
velocity held into incompressible and compressible components, that is, u = + u c , where 

V • = 0 and V x u c = 0. The DNS data base is investigated to ascertain whether the 

dilatational component u c is explicitly responsible for the reduced production. The Helmholtz 
decomposition leads to the following expression for the Reynolds shear stress anisotropy, 


612 = (1 - Xk)b{ 2 + Xkb^ (22) 

where Xk = K c / K is the compressible fraction of turbulent kinetic energy, 6[ 2 = u\' U 2 /2K 1 , 
and 6^2 = (m(m 2 — u{' U 2 )/2K c . It is found that the last term on the right hand side of Eq. (??) 
is much smaller than the first term on the right hand side primarily because Xk < 10% in all 
the simulations considered here. This implies that the reduced 612 of the solenoidal velocity 
component is responsible for the reduction in the overall value of 612 at high gradient Mach 
numbers. 

Thus, the DNS results of series A indicate that the large reduction of the turbulence growth 
rate that occurs when the gradient Mach number M g 0 increases is primarily due to the reduced 
Reynolds shear stress anisotropy of the solenoidal velocity component and not due to the explicit 
dilatational terms. 

The gradient Mach number is changed in series A by changing the shear rate S . A legitimate 
concern is that increasing S in the incompressible case may lead to decreased turbulence growth 
rates and the observed reduction of A in series B may not be a ‘compressibility effect’ after 
all. In order to address this concern, four cases AI1,AI2, AI3 and AI4 were performed where 
the incompressible Navier-Stokes equations were simulated. The initial data and parameters 
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were identical between All and Al, between AI2 and A2, between AI3 and A3, and between 
AI4 and A4 except that the incompressible Navier- Stokes equations were numerically solved in 
cases AI1-AI4 in contrast to the simulation of the compressible Navier- Stokes equations in cases 
A1-A4. Since the speed of sound c = oo in the incompressible case, the gradient Mach number 
and turbulent Mach number are both zero in cases AI1-AI4. The initial values of SK/e in cases 
All, AI2, AI3 and AI4 are 1.8, 3.6, 5.4, and 10.8 respectively. Fig. 5 shows the evolution of 
turbulent energy growth rate A in the incompressible simulations. The simulations had to be 
terminated by St = 12 so as to preserve good resolution of the small scales. The difference 
between the long-time growth rates among cases AI1-AI3 is small and, furthermore, the higher 
shear case AI3 seems to have a slightly larger value of A at the end of the simulation than the 
low shear case All contrary to the dramatic reduction with increased shear observed in the 
compressible simulations. However, the highest shear case AI4 has a slightly reduced growth 
rate probably because the initial value of SK/e = 10.8 is large enough to cause some rapid 
distortion effects in the turbulence evolution. Overall, it is clear that the reduced turbulence 
growth rate seen in series A of the simulations is indeed a compressibility effect. 

4.2 DNS results on the turbulence growth in series B 

The effect of increasing M t o, keeping M g q constant, on the turbulence growth rate is inves- 
tigated in series B of the simulations. According to the DNS, A evolves to asymptotic values 
that show a systematic decrease with increasing M t o- Fig. 6 shows that the high M t o case B3 
has an asymptotic value A^ ~ 0.11 which is smaller than the corresponding value A^ ~ 0.15 
in the low M t o case Bl. Although significant, the reduction is not as large as in series A where 
the high M g o case A4 has a growth rate A^ = 0.03 which is only 1/3 of that in case Al. 

Equation (??) is again used to distinguish between the contribution of production, dissipa- 
tion and explicit dilatational terms to the observed ‘compressibility effect’ of reduced turbulence 
growth rate. Figs. 7 and 8 show the evolution of normalized production — 2&i2, and the relative 
dissipation e s /V , respectively. Comparison of cases Bl and B3 in Figs. 7 and 8 show that case 
B3 has a lower production level and higher dissipation level relative to case Bl. This difference 
leads to the lower value of A for case B3 relative to case Bl in Fig. 6. Figure 9 shows the 
evolution of the pressure-dilatation and compressible dissipation. According to Fig. 9, the cases 
with higher M t o have larger dilatational contributions to the normalized growth rate. However, 
the difference in the dilatational terms between the various cases in Fig. 9 is small. 
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Thus, the compressibility effect of decreased growth rate of turbulent kinetic energy in series 
B is due to both a reduced level of production and an increase in the relative contribution of 
the turbulent dissipation rate. This is in contrast with series A of the simulations where the 
reduced level of production is solely responsible for the reduced growth rate of turbulent kinetic 
energy. 

4.3 Consolidated results from series A and B 

The growth rate of the turbulent kinetic energy ranges from a high of 0.15 in case B1 to 
a low of 0.03 in case A4. Thus, compressibility has a strong ‘stabilizing’ effect on the turbulence 
growth. This occurs in spite of rather small values of the dilatational component of velocity; 
the compressible fraction of turbulent kinetic energy \K is always less than 10% at the end of 
each simulation. However, the influence of compressibility on the pressure held is much more 
significant as shown below. A method to quantify compressibility effects on the pressure is 
provided by the pressure decomposition (Erlebacher et. al. 1990, Sarkar 1992) which splits the 
pressure held at a given time into p' = p 1 ' + p G \ where the incompressible pressure p 1 ' satishes 

VV' = -2 U hJ ('pu I /) g - 2U it i(pu I /) t j - (puf'u 1 /)^ - 2 Uijjipu 1 /) (23) 

while the remainder p G ' is the compressible pressure. Thus, Eq.(??), which is a Poisson equation 
for p l \ has right hand side terms that involve p and is free of p' - terms. Although the 
pressure decomposition is not unique unlike the corresponding Helmhlotz decomposition for 
the velocity, it is a useful method to quantify deviations of the pressure held from the baseline 
incompressible behavior. In the simulations, since p' is available from the compressible Navier- 
Stokes solution and p 1 ' is evaluated from Eq. (??), we obtain p G ' as the difference p' — p 1 ' . Figs. 
10 and 11 show the evolution of the ratio p^ms/Prins i n series A and B, respectively. It is clear 
that the ‘compressible pressure’ huctuation p unlike the dilatational velocity huctuation, is 
comparable to and actually larger than the incompressible counterpart. Figs. 12 and 13 show 
the evolution of the correlation coefficient R(pc,pi) = p G> p 1 ' / (PrmsPrms) ■ A systematic increase 
in the magnitude of R(pc,pi ) from the low M g q case Al to the high M g q case A4 is evident in 
Fig. 12. Also, the magnitude of R(pc,pi ) is generally larger in series A relative to series B. It 
appears from Figs. 10-13 that compressibility has a signihcant quantitative effect on the pressure 
held. The consequent effect on the pressure-strain correlation in the Reynolds stress transport 
equations could lead to the reduction in the magnitude of the anisotropy of the Reynolds shear 
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stress & 12 , and the observed increase (not discussed in detail here) in the magnitude of the 
normal stress anisotropy. Thus, pressure-strain models that are used for incompressible flows 
may not be applicable to high-speed shear flows. 

The gradient Mach number and turbulent Mach number evolve in time because of the 
temporal evolution in the transverse length scale / and the kinetic energy K, respectively. A 
plausible locality hypothesis for the influence of Mach number is that, after the initial transient 
and for sufficiently high Reynolds number, the value of either M g or M t or both at a given 
time determines the value of turbulence growth rate A at that time. From Figs. 1 and 6, it 
appears that A eventually becomes approximately constant in time. Therefore, if either M g or 
M t show a similar tendency to become invariant with time in the DNS, the locality hypothesis 
about the influence of Mach number would be supported. Figs. 14 and 15 show the evolution 
of M t in series A and B, respectively. Of the seven cases, only case A4 shows a possible 
tendency of M t to become constant. M t exhibits a monotone increase with time in all other 
cases. Thus, the DNS data on the long-time behavior of M t do not support the hypothesis 
that the instantaneous value of M t determines the turbulence growth rate. Figs. 16 and 17 
show the evolution of M g in series A and B, respectively. According to Fig. 16, the M g (St) 
curve tends to flatten with increasing St in cases A1-A3. Fig. 17 shows that after an initial 
increase in the series B simulations, M g decreases slightly for large time. Since the large eddies 
are constrained by the finite dimensions of the computational box, / stops increasing eventually 
and M g = SI /c stops increasing. The DNS data of Figs. 16 and 17 show a tendency for M g to 
become eventually constant and provide some support for the hypothesis that the turbulence 
growth rate depends on the instantaneous value of the gradient Mach number M g . Assuming 
that such a hypothesis is meaningful, we plot A^ normalized by the incompressible value of 
0.15 in our DNS against the terminal value of M g in Fig. 18. According to Fig. 18, there is 
a large reduction of the normalized turbulence growth rate in homogeneous shear flow when 
M g increases. The similarity of the curve in Fig. 18 with the variation (Kline, Cantwell and 
Lilley 1982, Papamoschou and Roshko 1988) of thickness growth rate db/dx with convective 
Mach number M c is evident. This similarity between the compressibility effect on the growth 
rate of turbulent kinetic energy in homogeneous shear flow and the compressibility effect on the 
thickness growth rate of the turbulent mixing layer is not surprising in view of Eq. (??) that 
relates the thickness growth rate db/dx to the kinetic energy growth rate A. 
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5. Implications of the DNS results for compressible shear flows 

The reduced growth rate of turbulent kinetic energy (and all its components) observed in 
the DNS of compressible homogeneous shear flow is analogous to the reduction in turbulence 
intensities (Elliott and Samimy 1990, Barre, Quine and Dussauge 1994) and thickness growth 
rate in the high-speed mixing layer. Inspection of Eq. (??) shows that the streamwise growth 
rate d8/dx of the mixing layer is proportional to A 00 / 6*2 where is the asymptotic growth 
rate of turbulent energy in the corresponding homogeneous shear flow problem, and C 2 defined 
by Eq. (??). Thus, the mixing layer growth rate d8/dx will decrease with increasing Mach 
number if either of the following factors occur, 

1. The asymptotic growth rate A^ of K in the corresponding homogeneous shear problem 
reduces. 

2. The coefficient C 2 increases, that is, the convection velocity U c becomes a relatively larger 
fraction of (U\ + tQ). 

According to our DNS, the first of the above conditions is met. Experimental or simulation 
data are required to investigate whether the second factor occurs. Thus, the decrease of A^ 
with the increase in Mach number that is observed in our DNS is consistent with the experi- 
mentally observed reduction in the growth rate of mixing layer thickness. More importantly, 
the phenomenon of reduced ‘efficiency’ of turbulence production by shear that occurs in the 
homogeneous shear flow DNS may be an important contributor to the reduced shear layer 
growth observed in experiments on the compressible mixing layer. A recent numerical study 
by Papamoschou and Lele (1993) of the evolution of the disturbance held of a small, isolated 
vortex placed in a hyperbolic-tangent mean flow profile has found that the growth rate of the 
perturbation energy decreases due to reduced production at high convective Mach numbers. 
Although the conditions (such as a small-amplitude disturbance and isolated vortex) consid- 
ered by Papamoschou and Lele (1993) are very different from those in the present DNS of 
homogeneous shear flow, it is interesting to note that both studies find that compressibility 
reduces the Reynolds shear stress. 

The effect of gradient Mach number M g appropriately defined using the mean compression 
rate D = | V • U\ instead of the shear rate S has been investigated by Jacquin, Cambon and Blin 
(1993) for the shock/turbulence interaction problem and by Cambon, Coleman and Mansour 
(1993) for the case of homogeneous compression in the rapid distortion theory (RDT) limit, 


16 



that is, for DK / e >> 1. Both investigations found that the growth of turbulent kinetic energy 
increases monotonically with the gradient Mach number, in contrast to the finding of the present 
paper. Of course, the present work differs in two crucial aspects from that of Jacquin et. al. 
(1993) and Cambon et. al. (1993) which may account for the different results. First, shear flow 
is considered here instead of compressed flow; and second, the mean shear rate is moderate so 
that non-linear effects are important in the present work while non-linear effects are negligible 
with respect to the large mean strain in the RDT approach used by the other authors. 

It is well known (Bradshaw 1977, Lele 1993) that compressibility effects over and above the 
influence of mean density are much stronger in the mixing layer relative to the boundary layer. 
The mean Mach number, either free-stream M ^ in the boundary layer or convective M c in 
the mixing layer, may take the same value, say 2; however, at this value of Mach number, the 
boundary layer does not show the large reduction in thickness growth rate that the shear layer 
does. This suggests that there may be another non-dimensional parameter based on the speed 
of sound which could conceivably distinguish between the boundary layer and the shear layer 
by taking different values although the mean Mach number is the same in these flows. We now 
show that the gradient Mach number M g could be such a parameter by estimating its value in 
the boundary layer and the mixing layer. 

Consider the adiabatic, zero pressure gradient boundary layer on a flat plate. The gradient 
Mach number, 


M g = 


( dll / dy) l 


(24) 


in the log-region, which occupies a large portion of the boundary layer, is estimated as follows. 
The Van Driest scaling for the law of the wall (Bradshaw 1977) in the adiabatic compressible 
boundary layer, 

rlTI . /7 Tt7 

(25) 


dU _ \/r w /p 
dy 


Ky 

where t w is the waU shear stress and k the von Karman constant, is used for the mean shear rate. 
Since the active, shear-stress producing motion scales with distance from the wall (Townsend 
1980), the approximate relation l — y is used for the integral length scale. Thus, Eq. (??) 
becomes 


M g 


Vt w /p 

KC 

M r 

K 


(26) 


17 



Here M T = \/t w /p w /c w is the friction Mach number based on wall friction and thermodynamic 


properties at the wall. The relation M r = M 00 ^C'//2 which links the friction Mach number to 
the free-stream Mach number is now substituted into Eq. (??) to obtain 


Mg ~ 


Mooy/Cj ] 2 


(27) 


Thus, the gradient Mach number is a constant given by Eq. (??) in the log-region of the 
high-speed boundary layer. 

We now consider the mixing layer. The mean velocity U in the mixing layer between two 
streams with velocities U\ and U 2 , respectively, is assumed (Schlichting 1979) to behave like an 
error-function, 

U-U 2 l + er%) 

FTM “ — 2 — (28) 

where r] = 2 (y — y c )/S(x ) is the similarity coordinate, y c is the position of the half- velocity 
point, S(x) is the mixing layer thickness, and erf ( 7 /) = (2 Jq exp( — u 2 )du. Let us estimate 
the value of M g at r] = 0 where the velocity gradient is maximum. From Eq. (??), the velocity 
gradient at r] = 0 is 2{JJ\ — U 2 ) / The length scale / of the turbulence is approximated 
by / ~ 8 which is consistent with the measurements of Wygnanski and Fiedler (1970) in an 
incompressible mixing layer, and the speed of sound at the center of the mixing layer is taken 
to be 0.5(ci + c 2 ). Then, the following expression 


M g ~ 2.2.1/ 


(29) 


is obtained as an estimate for the gradient Mach number at the centerline. 

We now compare the variation of M g as a function of M ^ and M c , respectively, for the 
mixing layer and the boundary layer. The Van Driest II expression for Cf(M 00 ) which agrees 
well with experimental data (Bradshaw 1977) is used in Eq. (??) to obtain M g as a function 
of Mqq in the boundary layer. The momentum thickness Reynolds number for the boundary 
layer is chosen to be Reg = 10 4 . The solid curve in Fig. 19 displays the variation of M g at the 
mixing layer centerline given by Eq. (??) while the dashed curve shows the variation of M g 
in the log-region of the boundary layer represented by Eq.(??). It is clear that the gradient 
Mach number M g for the mixing layer becomes much larger than that for the boundary layer 
when the mean Mach number increases. The smallness of M g in the boundary layer relative to 
the shear layer is a crucial difference between the two flows that could account for the much 
smaller compressibility effect in the boundary layer relative to the shear layer. The turbulence 
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production level is observed to decrease in our DNS when M g increases. Thus, on the basis of 
these considerations, the compressibility effect that reduces turbulence production is expected 
to be far smaller in a boundary layer relative to the mixing layer even if the mean flow Mach 
numbers have the same value in the two cases. This expectation is borne out by experimental 
results for the two flows. 

6. Conclusions 

DNS of a compressible fluid in turbulent homogeneous shear flow is performed for several 
cases in order to investigate and clarify some aspects of the influence of compressibility on 
turbulence. The two Mach numbers relevant to homogeneous shear flow - the turbulent Mach 
number M t and the gradient Mach number M g - were prescribed to have different initial values 
in the simulations. The initial value of M g is varied keeping the initial value of M t constant in 
series A of the simulations, and vice-versa in series B. Although the Reynolds number Re\ ~ 40 
at the end of the simulations is moderate, it is large enough for the flow to be turbulent. A total 
of six 128 3 simulations are performed. The simulations are used to investigate the influence 
of compressibility on the normalized growth rate of the turbulent kinetic energy K given by 
A = (1 / S K)(dK / dt) where S is the mean shear rate. The behavior of A is important because 
it is shown here through an analogy to be proportional to the thickness growth rate of the 
corresponding inhomogeneous mixing layer. The DNS results show that homogeneous shear 
flow, which is perhaps the simplest example of a turbulent shear flow, is strongly affected by 
compressibility. The asymptotic growth rate A^ of the turbulent kinetic energy decreases 
when either the initial gradient Mach number M g o or the initial turbulent Mach number M t o is 
increased, and furthermore, the reduction is substantially larger in the situation where M g o is 
increased. This ‘stabilizing’ effect of compressibility is analogous to the reduction in turbulent 
intensities and thickness growth rate observed in experiments on the high-speed mixing layer. 

The asymptotic growth rate of the turbulent kinetic energy depends on the eventual balance 
between the turbulent production by shear, the turbulent dissipation, and dilatational effects. 
A systematic comparison of these terms in the different cases is performed in order to identify 
probable causes for the reduced growth rate. Such a comparison leads to the conclusion that the 
reduction in the Reynolds shear stress anisotropy 612 and consequent reduction in the turbulence 
production level is predominantly responsible for the reduced growth rate of the turbulent 
kinetic energy. The dilatational terms (pressure dilatation and compressible dissipation) are 
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not negligible in the sense that they are observed to be as large as 20% of the turbulent kinetic 
energy budget. However, the difference in the values of the growth rate A between different 
cases is not due to differences in the dilatational terms. 

The gradient Mach number M g appears to be an important parameter in compressible shear 
flows. The turbulent energy growth rate decreases significantly when the gradient Mach number 
increases as shown in Fig. 18. For example, the growth rate is only half of its incompressible 
value at M g = 2. We have estimated the value of M g as a function of the convective Mach 
number M c in the mixing layer and as a function of the free stream Mach number M ^ in the 
wall boundary layer. The parameter M g increases significantly more rapidly in the mixing layer 
relative to the boundary layer, when the mean Mach number of the flow (M c in the mixing layer 
or Mqo in the boundary layer) increases. For example, a supersonic shear layer with M c = 1 
has M g ~ 2.2, while a supersonic boundary layer with M ^ = 2, Reg = 10 4 has M g ~ 0.16. It 
is well known that explicit compressibility effects over and above the effect of variable mean 
density are much larger in the supersonic shear layer relative to the supersonic boundary layer. 
Since the reduction in turbulence production by the mean shear and the consequent ‘stabilizing’ 
effect of compressibility in turbulent shear flows appears to be a strong function of the gradient 
Mach number in our DNS study, a possible reason for the difference between the extent of 
compressibility effects in the high-speed mixing layer and boundary layer is the large difference 
in values of the parameter M g . 
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